PERSISTENCE IN FLUCTUATING ENVIRONMENTS FOR 
INTERACTING STRUCTURED POPULATIONS 
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Abstract. Individuals within any species exhibit differences in size, develop- 
mental state, or spatial location. These differences coupled with environmental 
fluctuations in demographic rates can have subtle effects on population persis- 
tence and species coexistence. To understand these effects, we provide a general 
theory for coexistence of structured, interacting species living in a stochastic en- 
vironment. The theory is applicable to nonlinear, multi species matrix models 
with stochastically varying parameters. The theory relies on long-term growth 
rates of species corresponding to the dominant Lyapunov exponents of random 
matrix products. Our coexistence criterion requires that a convex combination 
of these long-term growth rates is positive with probability one whenever one 
or more species are at low density. When this condition holds, the community 
is stochastically persistent: the fraction of time that a species density goes be- 
low 5 > approaches zero as 5 approaches zero. Applications to predator-prey 
interactions in an autocorrelated environment, a stochastic LPA model, and spa- 
tial lottery models are provided. These applications demonstrate that positive 
autocorrelations in temporal fluctuations can disrupt predator-prey coexistence, 
fluctuations in log-fecundity can facilitate persistence in structured populations, 
and long-lived, relatively sedentary competing populations are likely to coexist in 
spatially and temporally heterogenous environments. 



1. Introduction 

All populations are structured and experience environmental fluctuations. Pop- 
ulation structure may arise to individual differences in age, size, and spatial loca- 
tion [Metz and Diekmann 1986 Caswell, 2001 , Holyoak et al. , 2005 . Temporal fluc- 
tuations in environmental factors such light, precipitation, and temperature occur 
in all natural marine, freshwater and terrestrial systems. Since these environmental 
factors can influence survival, growth, and reproduction, environmental fluctuations 
result in demographic fluctuations that may influence species persistence and the 
composition of ecological communities [Tuljapurkar 1990, Chesson, 2000b, Kuang 



and Chesson, 2009|. Here we present, for the first time, a general approach to 



studying coexistence of structured populations in fiuctuating environments 

For species interacting in an ecosystem, a fundamental question is what are the 
minimal conditions to ensure the long-term persistence of all species. Historically, 
theoretical ecologists characterize persistence by the existence of an asymptotic equi- 
librium in which the proportion of each population is strictly positive [May 1975 



Roughgarden 1979 



More recently, coexistence was equated with the existence 

a definition that 



Hastings, 1988 



of an attractor bounded away from extinction 
ensures populations will persist despite small, random perturbations of the popu- 
lations [Schreiber 2006, 2007] . However, "environmental perturbations are often 
vigourous shake-ups, rather than gentle stirrings" [Jansen and Sigmund 1998]. To 
account for large, but rare, perturbations, the concept of permanence, or uniform 



persistence, was introduced in late 1970s Freedman and Waltman, 1977 Schuster 
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et al. , 1979| . Uniform persistence requires that asymptotically species densities re- 
main uniformly bounded away from extinction. In addition, permanence requires 
that the system is dissipative i.e. asymptotically species densities remain uniformly 
bounded from above. Various mathematical approaches exist for verifying perma- 



nence [Hutson and Schmitt 1992 Smith and Thieme, 2011 including topological 
characterizations with respect to chain recurrence fButler and Waltman 1986 Hof- 



bauer and So 1989|, average Lyapunov functions [Hofbauer 1981 Hutson| 



Garay and Hofbauer 2003 , and measure theoretic approaches Schreiber 



1984 



2000 



Hofbauer and Schreiber, 20101. The latter two approaches involve the long-term, 



per-capita growth rates of species when rare. For discrete-time, unstructured mod- 
els of the form x^_,_^ = fi{xt)xl where xt = (xj, . . . ,x5^) is the vector of population 
densities at time t, the long-term growth rate of species i with initial community 
state xq = X equals 



1 

ri{x) = lim sup - Y^ log fi ( 



Garay and Hofbauer [2003| showed, under appropriate assumptions, that the system 
is permanence provided there exist weights pi, . . . ,pn associated with each species 
such that Y2iPi'''ii^) > for any initial condition x with one or more missing species 
(i.e. Hi ^* ~ 0)- Intuitively, the community persists if on average the community 
increases when rare. 

The permanence criterion for unstructured populations also extends to structured 
populations. However, in this case, the long-term growth rate is more complicated. 
Consider, for example, when both time and the structuring variables are discrete: 
xl_^_l = Ai{xt)xl where xl is a vector corresponding to the densities of the stages 
of species i, xt = (xj,...,x"), and Ai is a non- negative matrix. Then the long 
term growth rate rj(x) of species i corresponds to the dominant Lyapunov exponent 
associated with the matrices Ai{x) along the population trajectory: 

nix) = limsup-logll^j(xi) . . .Ai{xi)\\. 



t-^oo 



t 



At the extinction state x = 0, the long-term growth rate rj(0) simply corresponds 
to the In of the largest eigenvalue of ^j(O). For structured single-species models. 
Gushing [1998 , Kon et al. [2004 proved that ri(0) > implies permanence. For 



structured, continuous-time, multiple species models, rj(x) can be defined in an 
analogous manner to the discrete-time case using the fundamental matrix of the 
variational equation. Hofbauer and Schreiber 12010 showed, under appropriate 



assumptions, that ^iPiri{x) > for all x in the extinction set is sufficient for 
permanence. For discrete-time structured models, however, there exists no general 
proof of this fact. However, Salceanu and Smith 2009a|b 2010 used long-term 
growth rates rj(x) to establish persistence results for some stage structured, discrete- 
time models. When both time and the structuring variables are continuous, the 
models become infinite dimensional and may be be formulated as partial differential 
equations or functional differential equations. Much work has been done is this 



direction Hutson and Moran, 1987 Zhao and Hutson 



Magal et al. 2010 Xu and Zhao 2003 Jin and Zhao 



1994, Thieme, 2009 2011 



20091. In particular, for 



reaction diffusion equations, the long-term growth rates correspond to growth rates 
of semi-groups of linear operators and, YliPiri{x) > for all x in the extinction 
set also ensures permanence for these models [Hutson and Moran 1987, Zhao and 
Hutsonl 119941 ICantrell and Cosnerl |2003| . 
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Environmental stochasticity can be a potent force for disrupting population per- 
sistence yet maintaining biodiversity. Classical stochastic demography theory for 
stochastic matrix models xj+i = A{t)xt shows that temporally uncorrelated fluctu- 
ations in the projection matrices A(t) reduce the long-term growth rates of popu- 
lations when rare Tuljapurkar, 1990, Boyce et al. , 2006 . Hence, increases in the 



magnitude of these uncorrelated fluctuations can shift populations from persisting 
to asymptotic extinction. Under suitable conditions, the long-term growth rate for 
these models is given by the limit r = limt_^oo j In ||^(t) . . . ^(1)|| with probabil- 
ity one. When r > 0, the population grows exponentially with probability one for 
these density- independent models. When r < 0, the population declines exponen- 
tially with probability one. Hardin et al. 19881 and Benai'm and Schreiber 12009 



under different sets of assumptions, proved that these conclusions extend to mod- 
els with compensating density-dependence. However, instead of growing without 
bound when r > 0, the populations converge to a positive stationary distribution 
with probability one. These results, however, do not apply to models with over- 
compensating density-dependence or, more generally, non-monotonic responses of 
demography to density. 



Environmental stochasticity can promote diversity through the storage effect Ches ■ 



son and Warner 1981, Chesson 1982 in which asynchronous fluctuations of favor- 
able conditions can allow long-lived species competing for space to coexist. The 
theory for coexistence in stochastic environments has focused on stochastic differ- 
ence equations of the form xl_^_l = xlfi{xt,Ct+i) where ^i,i^2, ■ • • are a sequence of 



independent, identically distributed random variables (for a review see Schreiber 
201 1| ). [Schreiber et al.[ 2011| prove that coexistence, in a suitable sense, oc- 
curs provided that ^^Piri{x) > with probability for all x in the extinction set. 
Similar to the deterministic case, the long-term growth rate of species i equals 
ri{x) = limsup4__^oo j Yls=i ^^Sfii^s)- Here, stochastic coexistence implies that each 
species spends an arbitrarily small fraction of time near arbitrarily small densities. 
Our main result implies that the "community increases when rare" persistence 
criterion also applies to interacting, structured populations in stochastic environ- 
ments. Our model, assumptions, and a definition of stochastic persistence are pre- 
sented in Section 2. Except for a compactness assumption, our assumptions are 
quite minimal allowing for overcompensating density dependence and correlated en- 
vironmental fluctuations. Long-term growth rates for these models and our main 
theorem are stated in Section 3. We apply our results to stochastic models of 
predator-prey interactions, stage-structured beetle dynamics, and competition in 
spatial heterogenous environments. The stochastic models for predator-prey inter- 
actions are presented in Section 4 and examine to what extent "colored" environ- 
mental fluctuations facilitate predator-prey coexistence. In Section 5, we develop 
precise criteria for persistence and exclusion for structured single species models and 
apply these results to the classic stochastic model of larvae-pupae-adult dynamics 
of flour beetles Costantino et al. , 1995, Dennis et al. , 1995, Costantino et al. 1997 



Henson and Gushing 



1989 Gyllenberg et a. 



1997 



and Botsford, 2006, Schreiber 



and metapopulation dynamics Harrison and Quinn 



2010 



19961 |Metz and GyllenbeTJl [200T| [Roy et al.t|2005| [Hastings 



We show, contrary to initial expectations, 
that multiple noise with logarithmic means of zero can facilitate persistence. In 
Section 6, we examine spatial-explicit lottery models [Chesson 1985 2000a|b[ to 
illustrate how spatial and temporal heterogeneity, collectively, mediate coexistence 
for transitive and intransitive competitive communities. 
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2. Model and assumptions 

We study the dynamics of m interacting populations in a random environment. 
Each individual in population i can be in one of rii individual states such as their 
age, size, or location. Let XI = {Xl^, . . . , X^"') denote the row vector of populations 
abundances of individuals in different states for population i. XI lies in the non- 
negative cone M"'. The population state is the row vector Xt = {X^, . . . ,X™) that 
lies in the non-negative cone M" where n = X^j^i^j- To account for environment 
fluctuations, we consider a sequence of random variables, £,1,^2, ■ ■ ■ ,Ct, ■ ■ ■ where ^t 
represents the state of the environment at time t. 

To define the population dynamics, we consider projection matrices for each pop- 
ulation that depend on the population state and the environmental state. More 
precisely, for each i, let Ai{X,S,) be a non- negative, rii x m matrix whose j~k-th. 
entry corresponds to the contribution of individuals in state j to individuals in state 
k e.g. individuals transitioning from state j to state k or the mean number of off- 
spring in state k produced by individuals in state j. Using these projection matrices 
and the sequence of environmental states, the population dynamic of population i 
is given by 

^t+i = XlAi{S,t+i,Xt). 

where X| multiplies on the left hand side of Ai(^t+i,Xt) as it is a row vector. If we 
define A{^, X) to be the n x n block diagonal matrix diag(^i(.^, X), . . . , ^m(C) ^)); 
then the dynamics of the interacting populations are given by 

(1) Xt+i = XtA{^t+i,Xt). 

For these dynamics, we make the following assumptions: 

HI: ^1,^2) •• • is an ergodic stationary sequence in a compact Polish space E. 

H2: For each i, {£,, X) 1— ;■ Aj(^, X) is a continuous map into the space of rii x ni 
non- negative matrices. 

H3: For each population i, the matrix Ai has fixed sign structure corresponding 
to a primitive matrix. More precisely, for each i, there is a n^ x rii matrix Pi 
whose entries are zeros and ones such that all the entries of P/ are positive 
for some i > 1 and for all (X,^) the j-k-th entry of Ai{^,X) equals zero if 
and only if j-kth entry Pi equals zero. 

H4: There exists a compact set S C M" such that for all Xq G W^, Xt & S for 
all t sufficiently large. 

Our analysis focuses on whether the interacting populations tend, in an appro- 
priate stochastic sense, to be bounded away from extinction. Extinction of one or 
more population corresponds to the population state lying in the extinction set 

Sq = {x(^RI\YI\\x'\\=0} 

i 

where ||x*|| = X]^=i ^*"' corresponds to the i^ norm of x*. Given Xq = x, we define 
stochastic persistence in terms of the empirical measure 

(2) n^ = I E ^^^ 

s=l 

where 5y denotes a Dirac measure at y. These empirical measures are random 
measures describing the distribution of the observed population dynamics up to 
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time t. In particular, for any Borel set B C S, 

is the fraction of time that the populations spent time in the set B. For instance, if 
we define 

Srj = {x £ 5'|||x*|| < e for some i}, 

then Ilf{Srj) is the fraction of time that the total abundance of some population is 
less than e given Xq = x. 

Definition 2.1. For all rj > 0, define 

Sr^ = {x £ S : ||x*|| < r] for some i}. 
M is stochastically persistent if for all e > 0, there exists rj > such that 

UUS,) < e 
for t sufficiently large and all x £ S \ Sq. 

The set S^ corresponds to community states where one or more populations have a 
density less than rj. Therefore, stochastic persistence corresponds to all populations 
spending an arbitrarily small fraction of time at arbitrarily low densities. 

3. Results 

3.1. Long-term gro^vth rates and a persistence theorem. Understanding per- 
sistence often involves understanding what happens to each population when it is 
rare. To this end, we need to understand the propensity of the population to increase 
or decrease in the long term. Since 

Xi = XlAi{^uXo)Ai{^2,Xi) . ..Ai{it,Xt-i), 

one might be interested in the long-term "growth" of random product 

Ai{ii,X^)Ai{e,2,Xi)...A,{^uXt-i) 

of matrices as f — )• oo. One measurement of this long-term growth rate when Xq = x 
is the random variable 

ri{x) = limsup - log || ^^(^1, Xo)^i(6, ^i) • • • ^*(6, ^t-i)ll 

t— 5>00 t 

Population i is tending to show periods of increase when ri{x) > and asymp- 
totically decreasing when rj(x) < 0. The limsupj_j.oo instead of limt_s.oo in the 
definition of rj(x) is, in general, necessary. However, as we discuss in Section 3.2, 
the lim sup^_^oo can be replaced by limi_5.oo on sets of "full measure" . 

An expected, yet useful property of rj(x) is that ri{x) < with probability one 
whenever ||x*|| > 0. In words, whenever population i is present, its per-capita growth 
rate in the long-term is non-positive. This fact follows from XI being bounded 
above for t > 0. Furthermore, on the event of {limsup(_>oo ll-^'^tll > 0}) ^^ S^t that 
ri{x) = with probability one. In words, if population z's density infinitely often is 
bounded below some minimal density, then its long-term growth rate is zero as it is 
not tending to extinction and its densities are bounded from above. Both of these 
facts are proven in the Appendices. 

Our main result extends the persistence conditions discussed in the introduction 
to stochastic models of interacting, structured populations. Namely, if the commu- 
nity increases on average when rare, then the community persists. 
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Theorem 3.1. // there exist positive constants pi, ■ ■ ■ ,Pm such that 
(3) y^Piri{x) > with probability one 

i 

for all X £ So, then ([I]) is stochastically persistent. 



3.2. A refinement using invariant measures. The proof of Theorem 3.1 follows 
from a more general result that we now present. For this result, we shows that one 
need not verify the persistence condition [3] for all points x in the extinction set Sq. 
It suffices to verify the persistence condition for invariant measures of the process 
which concentrate on the extinction set. 

Definition 3.2. A Borel probability measure ^ on E x S an invariant measure for 
(II]) provided that 

(i) P[Ct eB]= fi{B X S) for all Borel sets B CE, and 

(ii) if F[{(,o,Xo) £ B] = pi{B) for all Borel sets B C E x S, then P[(6,^t) G 
B] = ix{B) for all Borel sets B C E x S andt>0. 

Condition (i) ensures that invariant measure is consistent with the environmental 
dynamics. Condition (ii) implies that if the system initially follows the distribution 
of n, then it follows this distribution for all time. When this occurs, we say (^t, X^) 
is stationary with respect to fi. One can think of invariant measures as the stochastic 
analog of equilibria for deterministic dynamical systems; if the population statistics 
initially follow /x, then it follows /i for all time. 

When an invariant measure ^ is statistically indecomposable, it is ergodic. More 
precisely, fi is ergodic if /U = a/ii + (1 — a)/i2 for < a < 1 and invariant measures 
/ij implies fii = /i2- 

Definition 3.3. // {(,t,Xt) is stationary with respect to jj,, the subadditive ergodic 
theorem implies that rj(Xo) is well-defined with probability one. Moreover, we call 
the expected value 

r,{f,) = jE[r,{Xo)\Xo = x,Ci = CMdx,dC) 

to be long-term growth rate of species i with respect to ^u. When /x is ergodic, 
the subadditive ergodic theorem implies that ri(Xo) equals ri{fj.) for ix-almost every 

(^0,6). 



With these definitions, we can rephrase Theorem 3.1 in terms of the invasion rates 
rj(//) as well as provide an alternative characterization of the coexistence condition. 

Theorem 3.4. If one of the following equivalent conditions hold 

(i) r=K(/u) := maxi<j<m?'j(Ai) > for every invariant probability measure with 

fi{So) = 1, or 
(ii) there exist positive constants pi, ■ . . ,Pm such that 

^Ptriin) > 

i 

for every ergodic probability measure with fi{So) = I, or 
(iii) there exist positive constants pi, ■ . . ,Pm such that 

2.Pi^i{x) > with probability one 

i 

for all x £ Sq 
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then (IT]) is stochastically persistent. 



To illustrate the applicability of Theorem 3.4, we apply the coexistence criterion to 
stochastic models of predator-prey interactions, stage-structured populations with 
over-compensating density-dependence, and transitive and intransitive competition 
in spatially heterogeneous environments. 



4. PrEDATOR-PREY DYNAMICS IN AUTO-CORRELATED ENVIRONMENTS 



For unstructured populations. Theorem 3.4 extends Schreiber et al. 2011 



s cri- 



teria for persistence to temporally correlated environments. These temporal corre- 
lations can have substantial consequences for coexistence as we illustrate now for a 
stochastic model of predator-prey interactions. In the absence of the predator, the 
prey, with density Nt at time t, exhibits a noisy Beverton-Holt dynamic given by 



(4) 



N, 



t+i 



Rt+iNt 
1 + aNt 



where Rt is a stationary, ergodic sequence of random variables corresponding to 
the intrinsic fitness of the prey at time t, and a > corresponds to the strength of 
intraspecific competition. To ensure the persistence of the prey in the absence of the 
predator, we assume that E[lni2i] > and E[lni?i] < oo. Under these assumptions. 
Theorem 1 of Benaim and Schreiber 2009] implies that Nt converges in distribution 



to a positive random variable N whenever Nq > 0. Moreover, the empirical measures 
nj ' with A^ > 0, P = converge almost surely to the law u of the random vector 
(A^, 0) i.e. the invariant probability measure satisfying i^{A) = P[(iV, 0) G A] for any 
Borel set AcM^. 

Let Pt be the density of predators at time t and exp(— 6P() be the fraction of 
prey that "escape" predation during generation t where b is the predator attack 
rate. The mean number of predators offspring produced per consumed prey is c, 
while s corresponds to the fraction of predators that survive to the next time step. 
Following May and Hassell, we get the following model of predator-prey dynamics 



(5) 



Nt+i 
Pt+i 



Rt+iNt 



l + aNt 
cNt{l — exp 



expi-bPt) 

hPt)) + sPt 



To ensures that (JHJ satisfies the assumptions of Theorem 3.1 we assume Rt takes 
values in the half open interval (0, i?*]. Since Nt+i < Rt+i/o- < R* /a and Pt+i < 
cNt + sPt < R* /a + sPt, Xt = {Nt,Pt) eventually enters and remains in the compact 
set 

S = [0, R*/a] X [0, eR*/{a{l - s))]. 



To apply Theorem 3.1, we need to evaluate ri({N,P)) for all A^ > 0,P > 

0. Since (0, Pt) converges to (0, 0) with probability one 
E[ln Rt] > and r2((0,P)) = Ins < 



with either A^ = or P 

whenever Pq > 0, we have ri((0,P)) 



whenever P > 0. Since IlJ ' with A^ > converges almost surely to u, we get 



ri((A^,0)) = by Proposition 8.18 and 



(6) 



r2((Ar,0))=E 



In cbN + s 



ln(c6x + s)u{dx) 
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By choosing pi = 1— e andp2 = e > for e sufficiently small (e.g. 0.5E[lni?i]/(E[lni?j] 
Ins)), we have Y2iPi'''iii-^j ^)) > ^ whenever NP = and 



(7) E In (cbN + s 



>0. 



Namely, the predator and prey coexist whenever the predator can invade the prey- 
only system. 

To see how temporal correlations influence whether the persistence criterion ([7]) 
holds or not, consider an environment that fluctuates randomly between good and 
bad years for the prey. On good years, Rt takes on the value rgood, while in bad 
years it takes on the value Rtad- Let the transitions between good and bad years 
be determined by a Markov chain where the probability of going from a bad year to 
a good year is p and the probability of going from a good year to a bad year is q. 
For simplicity, we assume that p = q in which case half of the years are good and 
half of the years are bad in the long run. Under these assumptions, the persistence 
criterion E[lniii] > for the prey is In {RgoodRbad) > 0. 

To estimate the left-hand side of (JTl), we consider the limiting cases of strongly 
negatively correlated environments (p ~ 1) and strongly positively correlated envi- 
ronments (p ~ 0). When p ~ 1, the environmental dynamics are nearly periodic 
switching nearly every other time step between good and bad years. Hence, one can 
approximate A^ by considering the positive, globally stable fixed point of 

-ti^good^t+l 

Xt+2 - ■— 

1 + axt+i 

_ RgoodRbadXt/jl + aXt)) 

1 + a{RbadXt/il + axt)) 

-^good-^badXt 

l + a{l + Rbad)xt 
which is given by ^n 1 r""^ \ • Hence, if p ~ 1, then A^ approximately puts half of 
its weight on ^"ni^'""'T^ and half of its weight on ^g°°f^'""'-^ and the persistence 
criterion ([7]) is approximately 



1 i__ / ,_ _RgoodRbad — 1 \ , ^ i / j, ^goodRbad — 1 



(8) - In fee -^;""-""" + g + 1 In fc^ -^"""-""" + g > q. 
^ ' 2 V a(l + Rbad) J 2 V «(1 + Rgood) J 

Next, consider the case that p ~ in which there are long runs of good years and 
long runs of bad years. Due to these long runs, one expects that half time N is near 
the value {Rgood — !)/« and half the time it is near the value nia,x{{Rhad — l)/0) 0}. 
If Rbad > 1; then the persistence criterion is approximately 

(9) I In f 6c^^^^^^^^ + s)+lln (bc^^^^^^^^ + s) > 
z\ a J 2 \ a J 

Relatively straightforward algebraic manipulations (e.g. exponentiating the left 
hand sides of ([S]) and ^ and multiplying by (1 -|- Ri,ad){^ + Rgood)) show that 
the left hand side of (Is]) is always greater than the left hand side of ([9]). 

Biological Interpretation 4.1. Positive autocorrelations, by increasing variabil- 
ity in prey density, hinders predator establishment and, thereby, coexistence of the 
predator and prey. In contrast, negative auto- correlations by reducing variability in 
prey density can facilitate predator-prey coexistence (Fig. Uy. 
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Figure 1 . Effect of temporal autocorrelations on predator-prey coexistence in a 

Markovian environment. In (a), the long-term growth rate r2((A'^, 0)) with 

A'^ > of the predator when rare is plotted as a function function of the 

temporal autocorrelation between good and bad reproductive years for the prey. 

In (b) and (c), the mean and interquartile ranges of long-term distribution of 

prey and predator densities are plotted as function of the temporal 
autocorrelation. Parameters: Rgood ~ 4, Rbad = 1.1, a = 0.01, c — 1, s — 0.1, 

6 = 0.01. 



5. Application to structured single species models 

For single species models with negative-density dependence, we can prove suf- 
ficient and necessary conditions for stochastic persistence. The following theorem 
implies that stochastic persistence occurs if the long-term growth rate ri(0) when 
rare is positive and asymptotic extinction occurs with probability one if this long- 
term growth rate is negative. 

Theorem 5.1. Assume that m = 1 (i.e. there is one species), H1-H4 hold and 
the entries of A(S^,x) = Ai{S,,x) are non-increasing functions of x. If ri{0) > 0, 
then 



(10) 



Xt+i — XtA{^t+i,Xt) 



is stochastically persistent. If ri{0) < 0, then limt-^oo Xt 
bility one. 



(0,0,..., 0) with proba- 



Proof. The first statement of this theorem follows from Theorem 3.1 



Assume that ri(0) < 0. Provided that Xq is nonnegative with at least one strictly 



positive entry, Ruelle's stochastic version of the Perron Frobenius Theorem iRuelle 



1979a, Proposition 3.2] and the entries of A{^,x) being non-increasing in x imply 

1 



hm llog\\Xt\\<l[m jlog\\XoA{CuO)...A{Ci,0)\\ 



with probability one. Hence, limi_s.oo Xt = (0, , 



r(0) < 
0) with probability one. 



D 



Theorem 5.1 extends Theorem 1 of Benaim and Schreiber [2009 as it allows 
for over-compensating density dependence and makes no assumptions about dif- 
ferentiability of X I— )• A{^^x). To illustrate its utility, we apply this result to the 
larvae-pupue-adult model of flour beetles and a metapopulation model. 



10 
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5.1. A stochastic Larvae-Pupae- Adult model for flour beatles. An impor- 
tant, empirically validated model in ecology is the "Larvae-Pupae- Adult" (LPA) 
model which describes flour beetle population dynamics Costantino et al. , 1995[ 



Dennis et al. 1995 Costantino et al. , 1997 . The model keeps track of the densities 



£t,Pt,at of larvae, pupae, and adults at time t. Adults produce b eggs each time 
step. These eggs are cannibalized by adults and larvae at rates Cea and c^i, respec- 
tively. The eggs escaping cannibalism become larvae. A fraction fii of larvae die at 
each time step. Larvae escaping mortality become pupae. Pupae are cannibalized 
by adults at a rate Cpa- Those individuals escaping cannibalism become adults. A 
fraction fXa of adults survive through a time step. These assumptions result in a 
system of three difference equations 



(11) 



£t+l = batexp{-Cel£t - Ceadt) 

Pt+l = (1 - fj-l)it 

at+i = {ptexp{-Cpaat) + (1 - /»„)«*) 



Environmental fluctuations has been included in these models in at least two 
ways. Dennis et al. 1995| assumed that each stage experienced random fluctuations 
due to multiplicative factors exp(^'),exp(^f),exp(^^) such that Q for i = l,p,a are 
normally distributed with mean zero i.e. on the log-scale the average effect of envi- 
ronmental fluctuations are accounted for by the deterministic model. Alternatively, 
Henson and Gushing [1997 considered periodic fluctuations in cannibalism rates 
due to fluctuations in the size Vt of the habitat i.e. the volume of the flour. In 
particular, they assumed that Cj = kj/\4 for i = l,a,p, for positive constants kj. If 
we include both of these stochastic effects into the deterministic model, we arrive 
at the following system of random difference equations. 



it+l = bat exp{-Kel£t/Vt+l - KgaOt/Vt+i + Cl) 

Pt+l = (1- w)^iexp(^f) 

at+i = (pt exp{-Kpaat/Vt+i) + (1 - fJ.a)at) exp(^") 



(12) 



Regarding these equations, we can use Theorem |3.4| to prove the following persistence 
result. When Q = with probability one for i = l,p,a, this theorem can be viewed 



as a stochastic extension of Theorem 4 of Henson and Gushing 1997 for periodic 

environments. 

Theorem 5.2. Assume c* > for i = l,p,a, fn G (0, 1) for i = l,a, ^f^t'^tJ '^^'^ 

Vt are ergodic and stationary sequences such that Q, log 14 £ {—M, M) for i = l,p,a 

and some M > 0, and (1 — fia) exp(^f ) G [0, 1 — 5] for some 5 > with probability 

one. Then there exists a critical birth rate bcrit > such that 

bcrit, then Xt = {it,pt,at) converges almost surely to 



Extinction: If b < 

(0,0,0) ast^oo. 
Stochastic persistence: If b > 

persistent. 

Moreover, if Q = ^f = ^f with probability one and¥,[Q 



bcrit, then the LPA model is stochastically 
0, then bcrit = /ia/(l-w)- 



Remark. The assumption that Q 



the normal distributions used by Dennis et al. 1995 



are compactly supported formerly excludes 
However, truncated normals 



with a very large M can approximate the normal distribution arbitrarily well. The 
assumption (1 — fia) exp(^^) S [0, 1 — 5] for some (5 > is more restrictive. However, 
from a biological standpoint, it is necessary as this term corresponds to the fraction 
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Figure 2. Effects fluctuations in fecundity and larval survival on the critical 
birth rate b required for persistence, ^j is normally distributed with mean and 
shown variance, £,t ~ £,t ~ ^ for all t and pa = 0.1034 (the value found in Table 
ID in 



Costantino at al 



1995 



of adults surviving to t he n ext time step. None the less, we conjecture that the 
conclusions of Theorem 5.2 hold when Q are normally distributed with mean 0. 

Theorem |5.2| implies that including multiplicative noise with log-mean zero has 
no effect on the deterministic persistence criterion when ^| = ^f = ^f with proba- 
bility one. However, when these random variables are not perfectly correlated, we 
conjecture that this form of multiplicative noise always decreases the critical birth 
rate (Fig. ^. To provide some mathematical evidence for this conjecture, we com- 
pute a small noise approximation for the invasion rate ri((5o) when population is 
rare iRuellel |1979al iTuljapurkarl 119901. Let 



Bt 



(l-w)exp(ef) 

^f 

Kb^l (l-;u„)exp(ef) 



be the linearization of the stochastic LPA model ( 12 ) at (L, P, A) = (0, 0, 0). Assume 
that Q = ^^ + eZl where ]E[Zj] = and IE[(Z^) ] = 1. A second order approximation 
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to Bf in e is 



Bt 



'0 (1-w) 

1 I (/ + ediag{eUf,Cf} + ^'diag{eUr,e?}V2 

.& (l-/i„) 



Let t" and w be the left and right Perron-eigenvectors of B such that ^- Vi = 1 and 
X^j fitfj = 1. Let vq be the associated Perron eigenvalue of B. Provided the Z^ are 
independent in time, a small noise approximation for the stochastic growth rate of 
the random products of Bt is 

~ r \ 2' 



(13) n(5o)«logro + ^ IE 



^ViWi{Ziy 



E 



y^ ■wj'Wi^t 



By Cauchy-Schwartz, 



E 



E 



ViWiZl 



< E 



< E 



E-^'(^)' 



Et;i'u;i(Zt*)' 



where the second line follows from < ViWi < 1 for all i. It follows that the order e^ 
correction term in ( 13 ) is non-negative and equals zero if and only if Zl = Zf = Z^ 



with probability one. Therefore, "small" multiplicative noise (with log-mean zero) 
which isn't perfectly correlated across the stages increases the stochastic growth 
rate and, therefore, decreases the critical birth rate bcru required for stochastic 
persistence. 

Biological Interpretation 5.3. For the LPA model, there is a critical mean fe- 
cundity, above which the population persists and below which the population goes 
asymptotically to extinction. Positive correlations between the survival rates of the 
different stages increases this critical mean fecundity. Negative correlations between 
the survival rates lower the critical mean fecundity. 

Proof of Theorem, 5.2. We begin by verifying H1-H4. HI and H2 follows from our 
assumptions. To verify H3, notice that the sign structure of nonlinear projection 



matrix At{^,X) for (12) is given by 



Since 



C 



C^ 






1 


0^ 








1 


1 





h 


^1 


1 


1 


1 


1 


2 


2 


1 


3 



At{^., X) has the sign structure of the primitive matrix C for all ^, X and t. Finally, 
to verify H4, define 



K = be 



2M-1 



/Ke 



Then 



^t+i < bat exp{-Keaat/Vt+i + ^l) < bat exp{-Keaat exp(-M) + M) < K 
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for alH > 0. Therefore, £j < i^ for t > 1. Hence, 



for all t > 2. Hence, 

at+i < Pte^ + (1 - S)at 

for all i > 2. Hence at < Ke /6 for t sufficiently large. The compact forward 
invariant set S = [0,K] x [0,Ke*^] x [0,Ke^^^]/5 satisfies H3. 
At low density we get 

/O (l-w)exp(er) 

Wt (l-/ijexp(en 

Define r(6) to be the dominant Lyapunov exponent of the random products of 
Bi,B2, ■ ■ . . Theorem 3.1 of Ruelle| [1979a implies that r{b) is differentiable for 
6 > and the derivative is given by (see, e.g., section 4.1 of Ruelle 1979a| ) equals 

vt{b)E3i'Wt+i{b) 



r'{b) = E 



>0 



vt{b)Bt{b)wt+i{b) 

where vt{b),wt{b) are the normalized left and right invariant sub-bundles associated 
with Bt{b) and E^i is the matrix with 1 in the 3—1 entry and entries otherwise. 
Since the numerator and denominators in the expectation are always positive, r{b) 
is a strictly increasing function of b. Since limfc_>o ^(^) = ~oo and limb_>oo ^{b) = oo, 
there exists bcrit > such that r(5) < for b < bcrit and r{b) > for b > bcrit 

If 6 > bcrit^ then r[b) > and Theorem 5.1 implies that (12) is stochastically 
persistent. On the other hand, if 6 < ftcrjtj then r{b) < and Theorem 5.1 implies 
that {it,pt,at) converges to (0,0,0) with probability one as t — )■ oo. 

The final assertion about the stochastic LPA model follows from observing that 
if ^f = ^1 = ^f with probability one for all t, then 

(1-w) 






Bt 



exp(^, 



(1 - ^^a) 



with probability one. Hence, r{b) = logro(6) + IE[q] where ?"o(6) is the dominant 
eigenvalue of the deterministic matrix 










1 

(l-^« 



Therefore, if E[log^|] = 0, then r{b) = logro(6). 
and Gushing [1997] showed that ro{b) > 1 if 6 > 



b < ^a/(l - w)- Hence, when ^J = ^f 
equals fj,a/{l — fJ-i) as claimed. 



Using Jury conditions 
/Xa/(1 - fii) and ro(6) 



Henson 



^^ with probability one and ^[Q 



< 1 if 

n 



5.2. Metapopulation dynamics. linteractions between movement and spatio- 
temporal heterogeneities determine how quickly a population grows or declines. 
Understanding the precise nature of these interactive effects is a central issue in 
population biology receiving increasing attention from theoretical, empirical, and 
applied perspectives [Petchey et al. 1997 Lundberg et al.', '2000', 'Gonzalez and Holt 

|2006j Hastings and Botsford 



2002 



2006 



Schmidt 2004 Roy et al. 2005 Boyce et al. 



Matthews and Gonzalez, 2007 Schreiber, 2010 . 
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A basic model accounting for these interactions involves a population living in 
a patchy environment with n patches. Let XI be the number of individuals in 
patch i at time t. Assuming Ricker density-dependent feedbacks at the patch scale, 
the fitness of an individual in patch i is Ajexp(— ajX^) at time t, where AJ is the 
maximal fitness and a, > measures the strength of infraspecific competition Let 
dij be the fraction of the population from patch i that disperse to patch j. Under 
these assumptions, the population dynamics are given by 

n 

(14) Xl^, = J2 dJ^KXl exp(-a,X,^) i = l,...,n. 



To write this model more compactly, let F{Xt,\t) be the diagonal matrix with 
diagonal entries Ai exp(— aiX^^), . . . , A^ exp(— OnX"), and D be the matrix whose 
i-j'th entry is given by dij. With this notation, (14) simplifies to 

Xt+i = XtF{Xt, \t)D 

Provided the A* are ergodic and stationary, \\ take values in a positive compact in- 
terval [A=K, A*] and D is a primitive matrix, then the hypotheses of Theorem |5 . 1| hold. 
In particular, stochastic persistence occurs only if ri(0), corresponding to the dom- 
inant Lyapunov exponent of the random matrix product F(0, \t)D . . . F{0, Xi)D, is 
positive. 

When populations are fully mixing (i.e. dij = 
derived a simple expression for ri(0) given by 



for all 



hJ) 



Metz et al. 



1983 



(15) 



ri(0) =E 



log E 



ViXl 



=1 



i.e. the temporal log-mean of the spatial arithmetic mean. Owing to the concavity 
of the log function, Jensens inequality applied to the spatial and temporal averages 



in (15) yields 



log iY,V^nK]] >]E 



log E 



ViXl 



> 



n 

E 



V^E[logXl 



The second inequality implies that dispersal can mediate persistence as ri(0) can 
be positive despite all local growth rates E[logAJ] being negative. However, the first 
inequality implies that dispersal-mediated persistence for well-mixed populations 
requires that the expected fitness IE[AJ] is greater than one in at least one patch. 
Schreiber [2010 developed first-order corrections to (15) for partially mixing popu- 



Vj lb e. A brief discussion of the implications of these results are 



lations where di 

provided in the Discussion. 

Biological Interpretation 5.4. Metapopulations with density-dependent growth 
can stochastically persist despite all local populations being extinction prone in the 
absence of immigration. 

6. Applications to competing species in space 
The roles of spatial and temporal heterogeneity in maintaining diversity is a funda- 



mental problem of practical and theoretical interest in population biology Chesson 
2000a|b[ [Loreau et al.[ |2003[ [Mouquet and Loreauj |2003[ [Davies et aTj |2005| . Tb 
examine the role of both forms of heterogeneity in maintaining diversity of compet- 
itive communities, we consider lottery-type models of m competing populations in 
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a landscape consisting of n patches. For there models, competition for vacant space 
determines the within patch dynamics, while dispersal between the patches cou- 
ples the local dynamics. After describing a general formulation of these models for 
an arbitrary number of species with potentially frequency-dependent interactions, 
we illustrate how to apply our results to case of two competing species and three 
competing species exhibiting an intransitive, rock-paper-scissor like dynamic. 

6.1. Formulation of the general model. To describe the general model, let X^"^ 
denote the fraction of patch r occupied by population i at time t. At each time step, 
a fraction e > of individuals die in each patch. The sites emptied by the dying 
individuals get randomly assigned to progeny in the patch. Birth rates within each 
patch are determined by local pair- wise interactions. Let ^^"'(r) be the "payoff" to 
strategy i interacting with strategy j in patch r at time t. Let 



t{r)=[cHr))_ 



l<i,j<m 

be the payoff matrix for patch r. The total number of progeny produced by an 
individual playing strategy i in patch r is YljCt -^V ■ Progeny disperse between 
patches with drs the fraction of progeny dispersing from patch s to patch r. Under 
these assumptions, the spatial-temporal dynamics of the competing populations are 
given by 

(16) Xti = , ^-^-^^-^^^'.^")-^^"^'^ + (1 _ e)xr-. 

Let Ai(^,X) be the matrix whose r — s entry is given by 

drsEjeHs)X^^ 



for r ^ s, and 



^E'sdrs'Ej,ia'Hs')X^''Xi^' 



for r = s. With these definitions, (16) is of the form of our models M. 

To illustrate the insights that can be gained from a persistence analysis of these 
models, we consider two special cases. The first case is a spatially explicit version of 
Chesson and Warner [1981 's lottery model. The second case is a spatial version of 



a stochastic rock-paper-scissor game. For both of these examples, we assume that a 
fraction d of all progeny disperse randomly to all patches and the remaining fraction 
1 — d do not disperse. Under this assumption, we get dgr = d/{m — 1) for s ^ r 
and dss = 1 — d. These populations are fully mixing when d = ^^^^^—— in which case 
d^r = — for all s, r. 



6.2. A spatially-explicit lottery model. The lottery model of Chesson and 



Warner 1981| assumes that the competing populations do not exhibit frequency 
dependent interactions. More specifically, the "payoffs" ^l'' (r) = Q{r) for all i,j are 
independent of the frequencies of the other species. Consequently, the model takes 
on a simpler form 



Esdrsms)X, 



(17) Xll^=e ^^-•'^^^.' ^ . +(i-g)xr. 

^ T.sdrsT.,et{s)Xr 
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where drs = i^^^zr^ for r ^ s and dgs = 1 — d. 

For two competing species (i.e. m = 2), define zi = (1, . . . , 1,0, . . . ,0) and 
Z2 = (0, . . . , 0, 1, . . . , 1) in which case the extinction set is 5o = {-zi, Z2}- Theo- 
rem [3]l] imphes that a sufficient condition for stochastic persistence is the existence 
of positive weights pi , p2 such that 

Piri{zi) +P2r2{zi) > and piri{z2) + P2r2{z2) > 

0. Hence, the persistence criterion 



r2[Z2) 



Proposition 8.18 imphes that ri{zi 
simphfies to 

ri{z2) > and r2{zi) > 0. 

In other words, persistence occurs if each species has a positive invasion rate when 
rare. 

To get some biofogical intuition from the mutual invisibility criterion, we con- 
sider the limiting cases of relatively sedentary populations (i.e. d ~ 0) and highly 
dispersive populations (i.e. d ^ 1). In these cases, we get explicit expressions for 
the realized per-capita growth rates ri{zj) that simplify further for short-lived (i.e. 
g: Ri 1) and long-lived (i.e. e w 0) species. Our analytical results are illustrated 
numerically in Fig. [3] 

6.2.1. Relatively sedentary populations. When populations are completely sedentary 
(i.e. d = ), the projection matrix A'^{^, (1, . . . , 1, 0, . . . , 0)) reduces to diagonal 
matrices whose r-th diagonal entry equals 

Mir) 



'^li 



+ 1 



The dominant Lyapunov exponent in this limiting case is given by 



r2{zij 



maxE 



log e 



'^Kr) 



Proposition 3 from Benaim and Schreiber [2009] implies that r2{zi) is a continu- 
ous function of d. Consequently, r2{xi) is positive for small d > provided that 



E 



log 



sM + i 



'(lir) 



is strictly positive for some patch r. Similarly, ri(z2) is pos- 



itive for small d > provided that E 



log 






is strictly positive for 



some patch r. Thus, coexistence for small d > occurs if 



maxE 

r 



lodeffl + l 



> and maxE 

r 



log £ 



iKr) 



+ 1 



>0. 



When e ~ 1 or e w 0, we get more explicit forms of this coexistence condition. 
When the populations are short-lived (e ~ 1), the coexistence condition simplifies to 
E[log a|(r)] > E[log a^(r)] and E[logai(s)] > E[logaj(s)] for some patches r ^ s. If 
the log- fecundities log a\{r) are normally distributed with mean ^i{r) and variance 
crf{r), this coexistence conditions becomes 

IJ-iir) > p.2{r) and /X2(s) > pi.i{s) 

for some patches r ^ s. In particular, coexistence requires that each species has at 
least one patch in which they have a higher geometric mean in their reproductive 
output. 
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When the populations are long lived (e ~ 0) and relatively sedentary (d ~ 0), the 
coexistence condition simplifies to 



E 



^Hr) 



> 1 and E 



(^ 



> 1 



for some patches r, s. Unlike short-lived populations, it is possible that both in- 
equalities are satisfied for the same patch. For example, when the log-fecundities 
log aj(r) are independent and normally distributed with mean /Ui(r) and variance 
(jf{r), the coexistence conditions become 



fi2{r) - fii{r) + 



> 1 



for some patch r and 



Mi(s) - ^J'2is) + 



atis 



+ a2(s) 



> 1 



for some patch s. Both conditions can be satisfied in the same patch r provided 
that (Ti(r) or a2{r) are sufficiently large. 

Biological Interpretation 6.1. For relatively sedentary populations, coexistence 
occurs if each species has a patch it can invade when rare. If the populations are 
also short-lived, coexistence requires that each species has a patch in which it is com- 
petitively dominant. Alternatively, if populations are long-lived, regional coexistence 
may occur if species coexist locally within a patch due to the storage effect. For un- 
correlated and log-normally distributed fecundities, this within-patch storage effect 
occurs if the difference in the mean log-fecundities is sufficiently smaller than the 
net variance in the log- fecundities. 

6.3. Well-mixed populations. For populations that are highly dispersive (i.e. 
d = ?^^^)j the spatially explicit Lottery model reduces to a spatially implicit model 
where 



riiS2 



E 



log e 



+ 1 



r2{Si) =E 
For short lived populations (e = 



log 

11, 






and 



these long-term growth rates simplify to 



ri{52 



r2{^i 



E 



E 






E 



E 






Since ri{d2) = —r2{5i), the persistence criterion that ri{62) > and r2{Si) > is 
not satisfied generically. 

Alternatively, for long-lived populations (e ~ 0), the invasion rates of well- mixed 
populations becomes to first-order in e > 0: 



r,{62)-e{E 
r2{6i)^e\E 
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1.010^ 



(I) 




Figure 3. Effects of dispersal and survival on coexistence of two species. The 

log-fecundities 5' are independent and normally distributed with means 

/xi = (5, 0, 5, 0, . . . , 0), /i2 = (0, 5, 0, . . . , 0) and variances cti = o-| = (1, . . . , 1) for 

(I) and (3, . . . , 3) for (II). The white lines correspond to the zero-lines of the 

respective Lyapunov exponents. 



We conjecture that this coexistence condition is less likely to be meet than the 
coexistence condition for relatively sedentary populations. To see why, we perform 
a small variance approximation of these invasion rates. Assume that ^l = ^^ +eZl{r) 
where Zl{r) are independent and identically distributed in i,r and E[Z^(r)] = for 
all i,r. Let a"^ = E[(Zj(r))^]. A second order Taylor's approximation in e yields 
the following approximation of the (rescaled) long-term growth rates for well- mixed 
populations 



(18) 



E 



Er&ir) 






and the following approximation for relatively sedentary populations 



(19) 



maxE 



CKr) 



1 






1 



Since (19) is greater than (18), persistence is more likely for relatively sedentary 



populations in this small noise limit. 

Biological Interpretation 6.2. Short-lived and highly dispersive competitors do 
not satisfy the coexistence condition. Long-lived and highly- dispersive competitors 
may coexistence. However, coexistence appears to be less likely than for sedentary 
populations as spatial averaging reduces the temporal variability experienced by both 
populations and, thereby, weakens the storage effect. 



6.4. The rock-paper-scissor game. In the last few years the rock-paper-scissor 
game, which might initially seem to be of purely theoretical interest, has emerged as 
playing an important role in describing the behavior of various real-world systems. 
These include the evolution of alternative male mating strategies in the side-blotched 



lizard Uta Stansburiana Sinervo and Lively, 1996] , the in vitro evolution of bac- 
terial populations [Kerr et al. 2002, Nahum et al.[ 2011| , the in vivo evolution of 
bacterial populations in mice [Kirkup and Riley 2004 , and the competition between 
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genotypes and species in plant communities [Lankau and Strauss 2007, Cameron 



et al. , 2009 . More generally, the rock-scissors-paper game - which is characterized 
by three strategies R, P and S, which satisfy the non-transitive relations: P beats R 
(in the absence of S), S beats P (in the absence of R), and R beats S (in the absence 
of P) - serves as a simple prototype for studying the dynamics of more complicated 



Leonard 1975 



Pascual, 2005 



non-transitive systems Buss and Jackson 1979, Paquin and Adams 1983, May and 



Schreiber, 1997, Schreiber and Rittenhouse, 2004, Vandermeer and 



Allesina and Levine 2011 . Here, we examine a simple spatial version 



of this evolutionary game in a fluctuating environment. 

Let xl{r),x^{r), and xf(r) be the frequencies of the rock, paper, and scissor 
strategies in patch r, respectively. All strategies in patch r receive a basal payoff of 
al at time t. Winners in an interaction in patch r receive a payoff of 6[ while losers 
pay a cost c^. Thus, the payoff matrix for the interacting populations in patch r is 



(r) 





bl 





bl 







We continue to assume that the fraction of progeny dispersing from patch r to patch 
s equals d/{m — 1) for s 7^ r and 1 — d otherwise. 

Our first result about the rock-paper-scissor model is that it exhibits a heteroclinic 
cycle in Sq between the three equilibria Ei = (1, . . . , 1,0, . . . ,0, 0, . . . , 0), E2 = 
(0, . . . , 0, 1, . . . , 1, 0, . . . , 0) and ^3 = (0, . . . , 0, 0, . . . , 0, 1, . . . , 1). 
Proposition 6.3. Assume d, e G (0,1] and al > cl, log a^ , log c^ , log 6^ G [—M,M] 
with probability one for some M > 0. // Xq > (0, . . . , 0) and Xq > (0, . . . , 0) and 
Xq = (0, . . . , 0), then limt_s.oo Xt = E2 with probability one. If Xq > (0, . . . , 0) and 
Xq > (0, . . . , 0) and Xq = (0, . . . , 0), then limt_s.oo xt = Ei with probability one. If 
Xq > (0, . . . , 0) and Xq > (0, . . . , 0) and Xq = (0, . . . , 0), then lim(_>oo xt = E3, with 
probability one. 

Proof. It suffices to prove the assertion for the case in which Xq > (0, . . . , 0) and Xq > 
(0, . . . , 0) and xg = (0, . . . , 0). Let 1 = (1, . . . , 1) G M". Our assumptions 6[ > and 
al > cl > imply there exists t/ > such that E'^{S,t+i, Xt) ^ exp{rj)E^ {^t+i, Xt) 
with probability one. It follows that 

1 1 

limsup - log \\Xt\\ = limsup - log \\E^{^t, ^t-i) • • • E^{^i,Xo)Xq\\ 



< limsup - log \\E\^t,Xt^i) . . . E^{^uXo)X'o\ 

t—>-oo I 

= limsup J log \\E\^t, Xt-i) . . . E\^,,Xo)Xi\ 



rj 
V 



< —rj 



where the last two lines follow from Proposition 8.14 and its Corollary. 



lim 



f— ^00 



XI 



as claimed. 



Hence, 

D 



Proposition 6.3 implies that all invariant measures on ^o are given by convex 



combinations of Dirac measures at the pure strategy equilibria. More specifically, 
define 6i to be the Dirac measure concentrated at the pure-strategy i equilibrium. 
Any invariant measure z^ on ^o is given by ^^Pi^i for a probability vector p = 
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{pi,P2,P3)- The persistence criterion of Theorem 3.4 requires that 



mm max 

peA3ie{i,2,3}' 



>0 



where A3 = {p £ M']_ '■ Pi + P2 + P3 = 1} the simplex of probabihty vectors. A 
standard algebraic calculation shows that this persistence criterion is satisfied if 
and only if 

r2{Si)r3{S2)rii63) > -r3((5i)ri(52)?"2(^3) 
i.e. the product of the positive invasion rates is greater than the absolute value of 
the product of the negative invasion rates. The symmetry of our model implies that 
all the positive invasion rates are equal and all the negative invasion rates are equal. 
Hence, coexistence requires 

'^3('^2) > -rii5i) 
As we did for the case of two competing species, we can derive more explicit 
coexistence criteria when the populations relatively sedentary (i.e. d ~ 0) or the 
populations are well-mixed (i.e. d = 1). For relatively sedentary populations, the 
coexistence criterion becomes 

^l + bl 



maxE 



log 1 -e + e- 



> — max E 



log 1 - e + e- 



For long-lived populations, this coexistence criterion simplifies further to 



maxE 



bl 



a 



> minE 



Alternatively, when the populations are well-mixed, the coexistence criterion be- 
comes 



E 



log 1 - e -F e 



Er«t 



> 



maxE 

r 



log 1 - e -F £ 



E.«[ 



/ JT 



For long-lived populations, this coexistence criterion simplifies further to 



E 



ErR 



> minE 

r 






Biological Interpretation 6.4. For relatively sedentary populations, coexistence 
only requires that average benefits (relative to the base payoff) in one patch is greater 
than the average costs (relative to the base payoff) in another patch. Negative cor- 
relations between benefits b^ and basal payoffs o[ promote coexistence. For highly 



dispersive species, If base payoffs are constant in space in time (i.e. 



a for 



all t, r), then coexistence requires the spatially and temporally averaged benefits of 
interactions exceed the spatially and temporally averaged costs of interactions. 

7. Discussion 

Understanding the conditions that ensure the long-term persistence of interacting 
populations is of fundamental theoretical and practical importance in population bi- 
ology. For deterministic models, coexistence naturally corresponds to an attr actor 
bounded away from extinction. Since populations often experience large perturba- 
tion, many authors have argued that the existence of a global attractor (i.e. perma- 
nence or uniform persistence) may be necessary for long-term persistence [Hofbauer 



and Sigmund 1998 Smith and Thieme, 2011 . Most populations experience sto- 



chastic fluctuations in their demographic parameters May, 1973 which raises the 
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question [May 1973, pg.621] "How are the various usages of the term [persistence] 
in deterministic and stochastic circumstances related?" Only recently has it been 
shown that the deterministic criteria for permanence extend naturally to criteria 



for stochastic persistence in stochastic difference and differential equations iBenaim 



et al. , 2008, Schreiber et al. 2011 . These criteria assume that the populations 
are unstructured (i.e. no differences among individuals) and environmental fluctu- 
ations are temporally uncorrelated. However, many populations are structured as 



highlighted in a recent special issue in Theoretical Population Biology Tuljapurkar 
et al. , 2012 devoted to this topic. Moreover, many environmental factors such 



as temperature and precipitation exhibit temporal autocorrelations IVasseur and 



Yodzis 2004 . Here, we prove that by using long-term growth rates when rare, 
the standard criteria for persistence extend to models of interacting populations 
experiencing correlated as well as uncorrelated environmental stochasticity, exhibit- 
ing within population structure, and any form of density-dependent feedbacks. To 
illustrate the utility of these criteria, we applied them to persistence of predator- 
prey interactions in auto-correlated environments, structured populations with over- 
compensating density-dependence, and competitors in spatially structured environ- 
ments. 



Mandelbrot 19821 proposed that environmental signals commonly found in na- 



ture may be composed of frequencies / that scale according to an inverse power 
law l/f^. With this scaling, uncorrelated (i.e. white) noise corresponds to /3 = 0, 
positively auto-correlated (i.e. red or brown) noise corresponds to /3 > 0, and nega- 
tively auto-correlated (e.g. blue) noise corresponds to /? < 0. Many environmental 
signals important to ecological processes including precipitation, mean air temper- 



ature, degree days, and seasonal indices exhibit positive (3 exponents Vasseur and 



Yodzis 20041. Consistent with prior work on models with compensating density 



dependence Roughgarden, 1975 Johst and Wissel, 1997, Petchey, 2000 , we found 



that positive autocorrelations in the maximal per-capita growth rate of species in- 
creases the long-term variability in their densities. If this species is the prey for a 
predatory species, we showed that this increased variability in prey densities reduced 
a predator's realized per-capita growth rate when rare. Hence, positive autocorrela- 
tions may impede predator-prey coexistence. In contrast, negative autocorrelations, 
possibly due to a biotic feedback between the prey species and its resources, may fa- 
cilitate coexistence by reducing variation in prey densities and, thereby, increase the 
predator's growth rate when rate. These results are qualitatively consistent with 
prior results that positive-autocorrelations in predator-prey systems can increase 



variation in prey and predator densities when they coexist Collie and Spencer 1994 



Ripa and Ives , 2003] . Specifically, in a simulation study of predator-prey interac- 



tions in pelagic fish stocks. Collie and Spencer 1994 found reddened noise resulted 



in predator-prey densities "to shift between high and low equilibrium levels" and, 
thereby, increase variability in their abundances. Similarly, using linear approxima- 
tions, Ripa and Ives 2003 found that environmental autocorrelations increased the 



amplitude of populations cycles. 

Classical stochastic demography theory [Tuljapurkar 1990, Boyce et al. , 2006| 
considers population growth rates in the absence of density-dependent feedbacks. 
Our results for populations experiencing negative-density dependence show that 
stochastic persistence depends on the population's long-term growth rate r(0) when 
rare. Hence, applying stochastic demography theory to r(0) provides insights into 
how environmental stochasticity interacts with population structure to determine 
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stochastic persistence. For example, a fundamental result from stochastic demogra- 
phy is that positive, within-year correlations between vital rates decreases r(0) and 
thereby may thwart stochastic persistence, a result consist with our analysis of the 
stochastic LPA model for flour beetle dynamics. Stochastic demography theory also 
highlights that temporal autocorrelations can have subtle effects on r(0). In particu- 
lar, for a density-independent version of the metapopulation model considered here, 



Schreiber| [2010] demonstrated that positive temporal autocorrelations can increase 
the metapopulation growth rate r(0) when rare for partially mixing populations, a 
prediction consistent with laboratory experiments [Matthews and Gonzalez 2007 



In contrast, Tuljapurkar and Haridas [2006 found that negative temporal autocorre- 



lations between years with and without fires increased the realized per-capita growth 

rate r(0) for models of the endangered herbaceous perennial Lomatium bradshawii. 

Spatial heterogeneity of populations has been shown theoretically and empiri- 



cally to have an effect on coexistence of competitive species (see e.g. Amarasekare 
[2003 or Chesson [2000b for a review). Coexistence requires species to exhibit 



Chesson 2000a 



niche differentiation that decrease the interspecific competition 
In a fluctuating environment, these niches can arise as differential responses to tem- 
poral variation |McGehee and Armstrong 1977, Armstrong and McGehee| 1980 



Chesson, 2000a|b[ , spatial variation May and Hassell, 1981, Chesson 2000a|b , or 



a combination of both forms of variation Chesson, 1985 . For the spatial lottery 
model where species disperse between a finite of patches and compete for micro sites 
within these patches, our coexistence criterion applies, and reduces to the mutual 
invisibility criterion. Although Chesson 1985 proved this result in the limit of an 



infinite number of patches with temporally uncorrelated fluctuations, our result is 
less restrictive as the number of patches can be small and temporal fluctuations can 
be autocorrelated. Using this mutual invisibility criterion, we derive explicit coexis- 
tence criteria for relatively sedentary populations and highly dispersive populations. 
In the former case, coexistence occurs if each species has a patch it can invade when 
rare. For short-lived populations, coexistence requires that each species has a patch 
in which it is competitively dominant. Alternatively, for long-lived populations, 
regional coexistence may occur if species coexist locally within a patch due to the 
storage effect Chesson and Warner, 1981, Chesson, 1982, 1994 in the one patch 



case. For highly dispersive populations, the coexistence criterion is only satisfied if 



populations exhibit overlapping generations, a conclusion consistent with Chesson 



1985 . By providing the first mathematical confirmation of the mutual invisibility 
criterion for the spatial lottery model with spatial and temporal variation, our re- 
sult opens the door for deeper investigations in understanding the relative roles of 
temporal variation, spatial heterogeneity, and dispersal on coexistence. Such inves- 
tigations could extend the work of Chesson [1985 and Muko and Iwasa 2000[ who 
showed that, without environmental variation, between-patch variation in the ratio 
of mortalities of two species promotes coexistence, while that of reproductive rates 
does not. 

For lottery models with three or more species, persistence criteria are more sub- 



tle and invasibility of all sub communities isn't always sufficient May and Leonard 



1975 



For example, in rock-paper-scissor communities where species 2 displaces 
species 1, 3 displaces 2 , and 1 displaces 3, all sub communities which consist of 



a single species are invasible by another, but coexistence may not occur [Hofbauer 



and Sigmund, 1998 Schreiber and Killingback, in press . For the deterministic 
models, coexistence requires that the geometric mean of the benefits of pair-wise 
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in 



interactions exceeds the costs of these interactions [Schreiber and Killingback 
press . Schreiber et al. [2011 and Schreiber and Kilhngback in press studied these 
interactions in models separately accounting for temporal fluctuations or spatial 
heterogeneity. In both cases, temporal heterogeneity or spatial heterogeneity can 
individually promote coexistence . Here we extend these result to intransitive com- 
munities experiencing both spatial heterogeneity and temporal fluctuations, thereby 
unifying this prior work. Our persistence criterion reduces to: the geometric mean 
of the positive long-term, low-density growth rates of each species (e.g. invasion 
rate of rock to scissor) is greater than the geometric mean of the absolute values 
of the negative, long-term, low-density growth rates (e.g. invasion rate of rock to 
paper). For relatively sedentary populations, coexistence only requires that aver- 
age benefits (relative to the base payoff) in one patch is greater than the average 
costs in another patch. Moreover, negative correlations between benefits and basal 
payoffs promote coexistence. For highly dispersive species, coexistence requires the 
spatially and temporally averaged benefits of interactions exceed the spatially and 
temporally averaged costs of interactions, assuming that base payoffs are constant 
in space and time. 

The theory of stochastic population dynamics is confronted with many, exciting 
challenges. While we have provided a sufficient condition for stochastic persistence, 
it is equally important to develop sufficient conditions for the asymptotic exclusion 
of one or more species with positive probability. In light of the deterministic theory, 
a natural conjecture in this direction is the following: if there exist non-negative 
weights pi, . . . ,pk such that 

^Pinifi) < 



for every ergodic measure fj, supported on the extinction set ^o, then there exist 
positive initial conditions such that Xt asymptotically approaches 5*0 with posi- 
tive probability. 



Benaim et al. 2008 



proved a stronger version of this conjecture 
for stochastic differential equation models where the diffusion term is small and 
the populations are unstructured. However, it is not clear whether there methods 
carry over to models with "large" noise or population structure. Another important 
challenge is relaxing the compactness assumption H4 for our stochastic persistence 
results. While this assumption is biologically realistic (i.e. populations always have 
an upper limit on their size), it is theoretically inconvenient as many natural models 
of environmental noise have non-compact distributions (e.g. log-normal or gamma 



distributions). One promising approach developed by Benaim and Schreiber 2009| 



for structured models of single species is identifying Lyapunov-like functions that de- 
crease on average when population densities get large. Finding sufficient conditions 
for "stochastic boundedness" is only half of the challenge, extending the stochastic 
persistent results to these "stochastically bounded" models will require additional 
innovations. Finally, and most importantly, there is a desperate need to develop 
more tools to analytically approximate or directly compute the long-term growth 
rates rj(^) when rare. 



8. Proof of Theorems 13.11 and 13.41 



Most of this Appendix is devoted to the proof of Theorem 3.4 from which The- 
orem 3.1 follows. In section 8.1, we recast our stochastic model (1) and our main 
hypothesis in Arnold's framework of random dynamical system Arnold, 1998 . Since 
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the dynamics we consider may not be invertible, in section 8.2 we extend the dy- 
namics to a dynamics on the set of possible trajectories to get an invertible dynamic 
on a larger space. Then we state Theorem 8.10 from which Theorem 3.4 follows. 



In section [8^ we prove basic results about the average per-capita growth rates Tj. 
In section |8.4[ we prove several basic results about occupational measures and their 



weak* limit points. A proof of Theorem 8.10 is provided in section 8.5 The proofs 



of Theorems 3.1 and 3.4 follows as a corollary of Theorem 8.10 (see Remark 8.11) 



8.1. Random dynamical systems frame^vork . To prove our main result, it is 
useful to embed (IT]) and assu mptions H1-H 4 within Arnold's general framework 
of random dynamical systems (Arnold 1998 ). Let Jl = E be the set of possible 
environmental trajectories, 3" = 6"^ be the product cr-algebra on ri, ^ : $7 i— )• il be the 
shift operator defined by 9{{ijjt}te'£) = {^t+i}tei^ arid Q be the probability measure 
on O satisfying 

Q({w G [7 : wt G ^0, • • ■,^t+k e Ek]) = P(Co G ^0, • • • ,efc e Ek) 

for any Borel sets Eq, . . . , En C E. Since £■ is a Polish space, the space Vt endowed 
with the product topology is Polish as well. Therefore, by the Kolmogorov con- 
sistency theorem, the probability measure Q is well defined, and by a theorem of 
Rokhlin Rokhlin 1964] , is ergodic with respect to Q. Randomness enters by choos- 
ing randomly a point uj = {ujt}ti^z G ^ with respect to the probability distribution 
Q and defining the environmental state at time t hy cot- 
In this framework, the dynamics (IT]) takes on the form 

/2QN f Xt+i{uj,x) = Xt{uj,x)A{uJt,Xt{u},x)) 

1 Xo{uj,x) 



xeS. 



We call (20), the random dynamical system determined by 



,A). 



To state the assumptions in this framework, define the skew product $ : x] 
0, X M'? associated with the dynamics (20) by 



and the projection maps pi : i^ x 
and p2{x,uj) = x. 



<^{uj,x) := {e{u}),xA{ujo,x)), 
— )• fi and p2 : i^ X 



i"' by pi{x,uj) = to 



Definition 8.1. A compact set K C 0, x M" is a global attractor for ^ if there 
exists a neighborhood V of K such that 

(i) for all {oj,x) G il x M", there exist T > such that ^^{uj,x) G V for all 

t>T; 
(ii) HV)CV andK = f]tm'^\y)- 
Assumptions H3-H4 take on the form 

H3': There exists a global attractor K C 0. xW^ ioi ^. 

H4': (7 is a compact space, Q is a Borel probability measure, and 9 is an invertible 

map that is ergodic with respect to Q, i.e. for all Borel set B C il., such that 

e-^{B) = B, we have Q{B) G {0, 1}. 

Assumptions H1-H2 do not need to be rewritten in the new framework. Assump- 
tion H3' is simply restatement of assumption H3 in the random dynamical systems 
framework. Since every ergodic stationary processes on a Polish space can be de- 
scribed as an ergodic measure preserving transformation (Kolmogorov consistency 
theorem and Rokhlin theorem), assumption H4' is less restrictive than H4. 
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To state Theorem 3.4 in this more general framework, we define invariant mea- 



sures for the random dynamical system (20). We foUow the definition given by 
First, recall some useful definitions and notations. We define int Ml := 



Arnold 



1998 



{x € Rl : Hi Xi > 0} and Mrf(M) the set of all d x d matrices over M. Let M be 
a metric space, and let T(M) be the space of Borel probability measures on M 
endowed with the weak* toplogy. If M' is also a metric space and f : M ^ M' is 
Borel measurable, then the induced linear map /* : 7{M) — ?• J'(M') associates with 
V e a'(M) the measure f*{v) G ^(M') defined by 

nv){B) = v{r\B)) 

for all Borel sets B in M' . li 6: M — )• M is a continuous map, a measure v G y(M) 
is called 6-invariant if z^(6'~^(i?)) = u{B) for aU Borel sets B £ M. A set i? C M is 
positively invariant if 0{B) C B. For every positively invariant compact set B, we 
call lia.v{6){B) the set of all ^-invariant measures supported on B. 

Definition 8.2. A probability measure fj, on 0, x Wl is invariant for the random 



dynamical system (20) if 

(i) ^ G Inv($), 

(ii) plifi) = Q, i.e. for all Borel sets D cQ, fiiD x W^) = 
For any positively invariant set il, x C where C C M" is compact, InV(Q(<I>)(ri x C) 
is the set of all measures /i satisfying (i) and (ii) such that ^{Vt x C) = 1. 

Remark 8.3. There is a bijection between the set of measures defined in Definition 



3.2 and the set InvQ($)(0 x M"). 



The following result is a consequence of Theorem 1.5.10 in Arnold 1998] . In fact, 
the topology defined in his definition 1.5.3 is finer than the weak* topology on the 
set of all probability measures on fi x C. 

Proposition 8.4. IfCc. W\. is a positively invariant compact set, then InvQ($)(rix 
C) is a nonempty, convex, compact subset ^[Vt x M"). 



The main assumption in Theorem 3.4 deals with the long-term growth rates which 



characterize, in some sense, the long term behavior of random matrix products (see 



Definition 3.3). In order to define those products in the new framework, we consider 
the maps Ai'.VlxS^ M„ . (M) defined by 

Ai{Lo,x) = Ai{ujo,x). 

While our choice of notation here differs slightly from the main text, this choice 
simplifies the proof. We write 

(21) A\{uj, x) := A,{uj, x)A^{^{u, x)) ■ ■ ■ Ai($*~^(a;, x)), 

with the convention that A^{ijj,x) = id, the identity matrix. 



Then, for each i G {1, . . . ^rn}, the asymptotic growth rate of the product (21) 
associated with (w, x) G il x R" is 

ri{x,ijj) := limsup - In ||A*(a;,x)||, 
t— >oo t 

which is finite due to assumption HI. According to Definition |8.2[ the invasion rate 
of species i with respect to an invariant measure fi G InV(Q(<I>) is 

ri{p) := I ri{uj,x)fi{du},dx). 

JnxR" 
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Given a point {uj,x) G il x M", let Ilti^,x) denote the empirical occupation 



measure of the process of ( 20 ) at time t defined by 

1 *"^ 
Ut{uj,x) := j^6x4co,x)- 

s=0 

For each Borel set B C M", the random variable Ilf{B) given by ([2|) is equal in 
distribution to the random variable Ilti-,x){B). 

For all T] > 0, recall that Srj := {x G W]^ : ||x*|| < rj for some i}. We can now 



rephrase Theorem 3.4 in the framework of random dynamical systems. 



Theorem 8.5. If one of the following equivalent conditions hold 

(i) r^:{fi) := maxo<j<m rj(^) > for every probability measure fj, E InvQ(<l>)(r2 x 

So), or 
(ii) there exist positive constants pi, ■ ■ ■ ,Pm such that 

^Pirii^i) > 

i 

for every ergodic probability measure /x E InvQ)(<I>)(il x Sq), or 
(iii) there exist positive constants pi, ■ ■ ■ ,Pm such that 

'^Piri{uj,x) > 

i 

for every x £ Sq and Q-almost all uj G Q, 
then for all e > 0, there exists r] > such that 

limsupnt(a;,x)(S',j) < e for Q- almost alluj, 

whenever x E M" \5'o. 

8.2. Trajectory space. The key element of the proof of Theorem|8.5|is Proposition 



8.13 due toiRuellei [1979bj in which it is crucial that dynamics induced by <1> is an 
homeomorphism. For now, it is not the case. Indeed, the map $ is, a priori, not 
invertible. The aim of this subsection is to extend the dynamics induced by ^ to 
an invertible dynamics on the larger set of possible trajectories. Then, we state an 



equivalent version of Theorem 8.5 in this larger space. 

By definition of the global attractor K, there exist a neighborhood V oip2{K) in 
M" such that $(17 x 1/) c Jl x F. By continuity of <1>, this inclusion still holds for 
the closure V of V , i.e. 

^{9. X F) Cl7 X F. 

This inclusion implies that, for every point (cj, x) E 12 x F, there exists a sequence 
{{ujt,xt)]ten C (ri X V)^ such that {ujq^xq) = {u:,x), and {u:t+i,xt+i) = $(wi,Xi) 
for all i > 0. Such a sequence is called a ^-positive trajectory. In order to create a 
past for all those ^-positive trajectories, let us pick a point {uj* ,x*) E x S\y , and 
consider the sequence space T := (0 x yu{x*})^ endowed with the product topology, 
and the shift operator : T -> T defined by Q{{{ujt,xt)]te-L) = {{'^t+i,xt+i)]tei- 
Since both and V U {x*} are compact, the space T is compact as well. 

Define 

r = {{{u:it,xt)}t(ii E T : {ujt+i,xt+i) = ^{ujt,xt), whenever xt E V}, 

endowed with the induced topology. Every ^-positive solution can be seen as an 
element of T by creating a fixed past (i.e. (wj, xt) = (w*, x*) for all t < 0). Then P is 
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the union of (i7 x {x*})^ and of the adherence of the set of all shifted (by 0* for some 
t G Z) <l'-positive solutions. Since F is a closed subset of the compact T, it is compact 
as well. Moreover F is invariant under 0, which implies that the restriction GL of 
O on F is well-defined. To simplify the presentation we still denote this restriction 
by ©. The projection map ttq : F — ;• $7 x 1/ U {x*} is defined by '/ro(7) = (wq, xq) for 
all 7 = {{ujt,xt)}t £ F. By definition, the map ttq is continuous and sujective. For 
now on, when we write 7 G F, we mean 7 = {(wt, Xt)}t&z- 
Next, we define the compact set of all ^-total trajectories as 

F+:=^o-i(J7xF), 

and the compact set of ^ -total- solution trajectory on the extinction set Sq as 

Fo:= ^0^^(17 x5o). 

The dynamic induced by <1> on il x F is linked to the dynamic induced by on F-|- 
by the following semi conjugacy 

(22) vro o e = $ o TTQ. 

For B <zT positively invariant and compact, define 

InvQ(e)(5) := {/i G Inv(e)(5) : p\ o ^*(/2) = Q}. 

Proposition 8.6. InvQ)(0)(F+) and InV(Q(0)(Fo) are compact and convex subsets 

o/?(r). 

Proof. Since F4. and Fq are positively invariant compacts, Inv(0)(F+) and InV(Q(0)(Fo) 
are non empty, compact and convex subsets of ^(F). Then, since Pi o ttq is contin- 
uous, InvQ(0)(F+) and InV(Q(0)(Fo) are compacts as closed subsets of y(F). The 
convexity of InV(Q(0)(F+) and InvQ(0)(Fo) is a consequence of the convexity of 
Inv(0)(F-t-) and Inv(0)(Fo), and the linearity of p^ o vtq. ■ 
As a consequence of equation ( 22 ) , we have 

Proposition 8.7. For every Q-invariant measure fx supported on F+, TTQ^fi) is $- 
invariant. 

Proof. Let /i be a 0-invariant measure supported on F+. Therefore the measure 
7rQ(/i) is supported by fi x F. Let B C il. x V he a Borel set. We have 

7:*{f,){^-\B)) = p{n,H<^-\B))) 

= fi{7:Q^{<^-^{B)r\nxV)) 

= f,{{nooe\^^)-\B)) 
= K^oHB)) 

The second equality is because the support of p, is indued in F_|_, and the third is a 
consequence of the conjugacy (22). ■ 

Corollary 8.8. 7rQ(InV(Q(0)(F_|-)) is a compact and convex subset ofl'n.VQ{^){il.xV). 



Proof. Since ttq is continuous and linear, Proposition 8.6 implies that 7rQ(InvQ)(0)(F_|_)) 
is a compact and convex. Proposition 8.7 implies that ■7r^InvQ)(0)(F+)) C InvQ)(<l>)(rix 

V). m 
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The map on F-i- can be seen as the extension of the map <I> on il x F. Now, we 
define the long-term growth rates for the extended dynamics induced by 0. 
For each i € {1, . . . , m}, define the maps Aj : F — )• M„. (M) by 

A-( ) = i M^*,x*) if xo = x* 
''^'^^ \ Ai{uJo,xo) either 

As ( |21[ ), we write 

(23) 4(7):=A.(7)---A.(0*-i(7)). 



The conjugacy (22) imphes that for all (w,x) € Q xV and all 7 G tTq {uj,x), we 
have 

(24) Ai{j) = Ai{co,x), 

for ah t > 0. 



Then the long-term growth rates for the product (24) is 

rj(7) :=limsup-lnp-(7)||, 
and, for a 0-invariant measure /2, the long-term growth rates is 



Proposition 8.9. For all species i, we have 

(i) ri{uj, x) = ri{'y), for all {to, x) ^ Q. x V and for all 7 G vTq" (cj, x), 
(ii) for all jx e InvQ(0)(F+), ■kUjI) € InvQ($)(r2 x V), and 

nijj) =ri{TT^{fl)). 



Proof. Assertion (i) is a consequence of equality ( 24 ) , and assertion (ii) is a conse- 



quence of assertion (i) and Corollary 8.8 



We can now state an equivalent version of Theorem 8.5 on the space of trajectories 
F. 

Theorem 8.10. If one of the following equivalent conditions hold 

(a) r^:{jl) := maxo<j<m rj(/i) > for every probability measure jl G InvQ(0)(Fo), 
or 

(b) there exist positive constants pi, ■ ■ ■ ,Pm such that 

^Piriifl) > 

i 

for every ergodic probability measure jl E InvQ(0)(Fo), or 

(c) there exist positive constants pi, ■ ■ ■ ,Pm such that 



'^Piri{uj,x) > 



for every x £ So and Q-almost all uj G 0,, 
then for all e > 0, there exists rj > such that 

limsupnt(a;,x)(S'^) < e for Q- almost alluj, 
whenever x £ M" \5'o. 
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Remark 8.11. Condition (c) of Theorem 8.10 and (Hi) Theorem 8.5 are equivalent, 
and the implications from conditions (Hi) to (ii) and (ii) to (i) of Theorem 8.5 are 
direct. The proof of Theorem 8.1(J\ (see section 8.5) shows that (a), (h) and (c) of 
Theorem 8.10 are equivalent. Finally, condition (i) of Theorem \8. 5| implies condition 
(a) of Theorem 8.10 as a direct consequence of assertion (ii) of Proposition \ 



Hence, Theorems 8.5 and 8.10 are equivalent. 



8.3. Properties of long-term gro'wth rates. In this section, we first state Propo- 



sition 3.2 of iRuelle [1979b| in its original framework, and extend it to ours. We use 
our extension to deduce some properties on tlie long-term growth rates which are 



crucial for the proof of Theorem 8.10 



Proposition 8.12 (Ruelle 1979b| ). Let E be a compact space, ^ : M+ xE ^ E be an 
homeomorphism, and vr : ExMr-^ E be the projection map defined by 7r(^, x) = ^. 
Consider a continuous map T : H — )• Mrf(M) and its transpose T* . Write 

T\() = T{0---T{^'~^0, 

and assume that 

A: for all ^ € E, r(0(M'^) C {0+} U intM^. 

Then there exist continuous maps u,v :E ^ M"^ with \\u{^)\\ = \\v{^)\\ = 1 such 
that 

(i) the line bundles E (resp. F) spanned by u{-) (resp. v{-)) are such that 

Rrf = S0 F-L where b G F{C)-^ if and only if (C, v{C)) = 0. 
(ii) E (resp. F) is T,"^ -invariant (resp. T*,^~^ -invariant), i.e. E{'i^{^)) = 

E{^)tI^) and F{^^)T*{'^^) = F{C), for all ^ £ E; 
(iii) there exist constants a < 1 and C > such that for all t > 0, and S, £ E, 

||6r(0 • • • T{¥-'0\\ < Ca'WaTiO ■ ■ ■ T{^'~'0\\, 
for all unit vectors a G E{^), b G F{^) . 

Since assumption H2 does not directly imply assumption A for the map Ai{-, •), 
we need to extend Ruelle's proposition to the case where 
Al': for ah ^ £ E, T{^) intM!^ C intM^^, and 

A2': there exists s > 1 such that, for all ^ G H, T{^) ■ • • r(*'*~i^)(M'^) C {0} U 
intM"^. 



8.12 



Proposition 8.13. The conclusions of Proposition 8.1S\ still hold under assump- 
tions Al'-A2'. 

Proof. Define the continuous map T' : S — )• H x Md(M) by 

r'(0 = T(0---r(vi/-i(e)). 

By assumption (A2'), T'{^)R'l C {0+}UintM^. Therefore, Proposition 
applies to the map T' and to the homeomorphism ^^ which give us maps u,v : 
M"^ with ||m(C)II = 11^(011 = 1) their respective vector bundles E{-),F{-), and some 
constants C,a verifying properties (i), (ii), and (iii). 

The vector bundles E{-),F{-) are our candidate bundles for T. We need only to 
check properties (ii) and (iii) for the map T as property (i) is immediate. 

We claim that 

xr*(*-*^) 



(25) 



lim ,, 



^011 



«(0, 
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uniformly on all compact of M'^ — {0+}. Before we prove (25), we show property 
(ii), i.e. £'(•) is T, ^-invariant, is a consequence (25). Let y S intM!^ — {0}, and 



^ G H. Continuity of T and equality (|25|) applied to y imply 

«(?)r(o 



lim ja^%ni) 



yr(vi/-*e)r*(^-*(^e)) 



lim 

t— >oo 



||yr*(*-*OII 

\\yT{^-'OT'{^-'m)\\ 



where the final line follows from ([25]) with C = ^C and x = yr('I'"*,^)/||yT(^"*^)|| 
which belongs to the compact {z £ M^ : ||2:|| = 1} for all t > 0. This proves property 
(ii) for E. The same argument for the transpose T'* implies property (ii) for F. Now 

— {0+} with ||a;|| = 1. For every t > 0, define st := i— [f]s 



we prove (25). Let x G S 

where [q] is the integer part of q. We have 

xT*(^-*0 = xT"*(^-*0r''^k^"*^"6- 
Since st < s for all t > 0, continuity of T(-), and assumption (Al') imply that there 
is a compact H C M^ — {0+} independent of x such that xr^*(^~*^) G H for all 



t > 0. Then, (25) is a consequence of inclusion (3.2) in the proof of Proposition 3.2 
in 



Ruelle 1979b| applied to the map T' . 



It remains to check property (iii): show that there exist a' ,C' > such that 
\\bT\^)\\ < C7V*||n(0r*(0ll for alH > s,^ G S,6 G F(0^. 
We have 

Since -F(-) is T*-invariant, bT'^^{^) G -F(^*'^)"'" and property (iii) for T' implies 



1 



|6r^*(0r'y(^'^'0ll < 



c{a-sY 



|n(e)r^nOr'y(^^'OII- 



The continuity of T{-,-) and u{-), and assumption Al' imply that there exist a 



constant R> such that 



max{||u;r'=(0|| : \\w\ 



-4<R. 



mm{l\u{i)T>'m : { e 2} 

for all A; < s and all ^ G H. Then property (iii) is verified with C" 



CR and 



a 



as . 



Proposition 8.13 applies to each continuous map Aj : F — )• M„,. (M), and to the 



homeomorphism Q on the compact space F. Then, for each of those maps, there 
exist row vector maps ni(-), Vi{-), their respective vector bundles Ei{-), Fj(-), and 



the constant Ci,ai > satisfying properties (i), (ii), and (iii) of Proposition 8.13 
For each i G {1, . . . , m}, define the continuous map Ci : F — )• M by 

C(7):=ln||u.(7)^.(7)ll- 



In the rest of this subsection, we deduce from Proposition 8.13 some crucial 
properties of the invasions rates. 

Proposition 8.14. For allj^T and every population i, ri{^) satisfies the following 
properties: 
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(i) 



ri(7) = limsup - In ||v^*(7)||, 
t— s>oo t 



(ii) 



for all V gM"'\{0}. 



ri 



1 *"^ 
(7)=limsup-J]C(e^(7)). 



s=0 



The proof is similar to that of proposition 1 in Hofbauer and Schreiber [2010 
Proof. Let 7 S F be fixed. To prove the first part, we start by showing that 

1 



(26) 



rj(7) =hmsup-ln||nj(7)^-(7)||. 



Let V G M"% u / 0. Since M"* = ^j 0F^(7), there exist a co nstan t a G M and a 
vector w G F^ such that v = aui{'-f) + w. Then, by Proposition 8.13, we have 

11^4(7)11 < a||n,(7)A*(7)ll + lkAK7)ll 

< \\u,{j)Ai{j)\\{a + aai\\w\\). 

Hence, 

hmsup - In ||u^*(7)|| < hmsup - In ||uj(7)74*(7)|| 

t— ^-oo f- t— >oo t 

for all w G M"'\{0}. Since ||A*(7)|| = sup||^||=;^ ||w^*(7)||, the last inequality implies 
that 

ri(7) < limsup - In ||'Uj(7)^-(7)|| < rii^j), 

which proves the equality (26). 



Now, we consider positive vector v G M^'\{0}. We show that the equality (26) is 
also satisfied for v. We write v = aui{'y) + w with a > and w G F,- (7). Proposition 



8.13 implies 



1^4(7)11 > a||n,(7)A*(7)||-||u;A*(7)ll 

> ||n,(7)A*(7)||(a-C,a*||u;||). 



Since a > 0, 



rti^) > limsup -In ||f A- (7) II > limsup - In ||iij(7)y4*(7)|| =ri{'y), 

t— >CXD t t— >00 t 

which completes the proof of assertion (i) . 

The second assertion results directly from the first assertion and the following 
equalities: 

ln||ni(7)^*+'(7)ll = In ||n.(7)A*(7)Ai(e*(7))|| 

= ln||n,(e(7)*)A,(G*(7))||||n,(7)4(7)|| 
= Ci(e*(7))+ln||n,(7)A*(7)||. 

The second step is a consequence of the invariance of the line bundle Ei . ■ 
RecaU that r+ = iTQ^i^ x V) and Tq = -n-Q^i^ x Sq). 

Remark 8.15. Assumption H2 implies that the sets Fq and F+\Fo are both posi- 
tively Q-invariant. Therefore every Q-invariant measure ji on F4. can he written as 
a convex combination of two Q-invariant measures i>o and v\ such that uo{Tq) = 1 
and i>i(F4_\Fo) = 1. 
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Corollary 8.16. For all 7 G r-|-\ro, and every i G {1, . . . ,m}, 

ni-f) < 0. 

Proof. Fix i £ {1, . . . ,m}, and 7 G r+\ro with {uj,x) := vro(7). By definition of 

r+\ro, x' G M"* and x* / 0. We have 

x*4(7) = x%(7)---Ai(e*-i7) 

= P2{<^\U},X)), 



where the second equality is a consequence of (23), and the third one follows from 
the definition of the cocycle $. Assumption H3' implies that there exists T > 
such that p2{^^{uj,x)) belongs to the compact set V for all t >T, which implies 
that there exists -R > such that ||x*A^(7)|| < R for all t > T. Assertion (i) of 



Proposition 8.14 applied to v = x* concludes the proof. ■ 

Now we give some properties of the invasion rate with respect to a G-invariant 
probability measure. 

Proposition 8.17. The invasion rate of each population i with respect to an Q- 
invariant measure jl satisfies the following property: 



Proof. This result is a direct consequence of property (ii) of Proposition 8.14 and 
the Birkhoff 's Ergodic Theorem applied to the continuous maps and C- ^ 

Proposition 8.18. Let jl he a Q-invariant measure. If fi is supported by r+\ro, 
then ri{fx) = for all i G {1, . . . , m}. 

Proof. Let jl be such a probability measure. Fix i G {1, . . . , ttt.}, and define the set 
p,r? ._ 1^ g Y_^ : ||p2(vro(0*(7)))*|| > ?/}• By assumption on the measure /2, there 
exists a real number 77* > such that /i(r*''^) > for all -q < rj* . 

The Poincare recurrence theorem applies to the map 0, and implies that for each 

T] < VI*, 

(27) /2({7 G r^'^'l e*(7) G F*'^ infinitely often }) = 1. 



Recall that the conjugacy (22) implies that for every 7 G F+ with 7ro(7) = (w,x), 
we have 

P2(^o(e*(7)))* = P2($*(vro(7)))' 
Then, equality (27) means that for /i-almost all 7 G F*''' with < rj < r]* , \\x'^Al{'y)\\ > 



rj infinitely often. Therefore, Proposition 8.14 (i), applied to f = x*, implies that 



rj(7) = for /i-almost all 7 G F*'^, with rj < r]* . Hence ri^-j) = for /i-almost all 
7 ^ Un>^ F*'-*^'" = F4.\Fo, which completes the proof. ■ 

— T)* 

8.4. Some Lemmas. Given a trajectory 7 G F+, the empirical occupation measure 
of{e*(7)}i>ois 

1 *"^ 

Ml) ■= j^^e^i-y)^ 

^ s=0 



STOCHASTIC PERSISTENCE OF INTERACTING STRUCTURED POPULATIONS 33 

and given a point {uj,x) £ Q xV, the empirical occupation measure of {^^{x,uj)}t 

t-i 



IS 



^ t-1 
^ s=0 



In this way, At{uj,x){0, x B) = Ilt{uj,x){B) for every Borel subset B C V, and 
xeV. 

Lemma 8.19. Let 7 G r+. Then for all t > we have 

7rS(At(7)) = At(vro(7)). 
Proof. Let B C Q x V he a Borel set, and 7 € F-i-. Then we have 

7r*(A,(7))(B) = Mj)i7r,HB)) 
1 *"^ 



^=0 

^ t-i 

At{Ml))iB). 



The third equahty is a consequence of the semi conjugacy (22) . ■ 

Lemma 8.20. There exists Q with Q{Cl) = 1 such that for all 7 G ttq" {Q x V), 
the set of all weak* limit point of the family of probability measures {At{'j)}t£N is a 
non-empty subset o/InvQ(0)(r+). 

Proof. Since Q is ergodic (assumption H4'), Birkhoff's Ergodic Theorem imphes 
that there exists a subset fi C il such that Q(J1) = 1, and for aU uj £ 0,, 

(28) }tJE'o 



Je^iui) 



s=0 



(in the weak* topology). Let (w, x) G J7 x F and 7 S vtq ^{u, x) C r+. For ah t G N, 
we have 



(29) 



1 *"^ 

plo^*CMi)) = ^^^ 



s=0 



Since r+ is positively 0-invariant and compact, the set of all weak* limit point of the 
family of probability measures {At(7)}tgN is a non-empty subset of J'(r+). Since the 
maps pi and ttq are continuous, equalities (28) and (29) imply that pi o 7rQ(/i) = Q. 
Moreover, Theorem 6.9 in Walters [1982| implies that fl is ©-invariant. Therefore, 
/i G InvQ(G)(r+), which concludes the proof. ■ 

Recall that Srj = {x £ S : ||x*|| < rj for some i}, and define the subset F^ := {7 G 
F+ : 7ro(7) G ^ x Sr,}. 



Lemma 8.21. // condtion (a) of Theorem 8.10 is satisfied, then for all e > there 
exists T]* > such that 

/i(F^)<e, V7/<r/*, 

for all jl G InvQ(0)(F+\Fo). 
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Proof. If false, there exist e > and a sequence of measures {/XnlneN C InvQ(0)(r+\ro) 



such that jjni^i/n) > £ foi' ^11 n > 1. By Proposition 8.6, let /i G InV(Q(0)(r+) 
be a weak* limit point of the sequence {/injneN- Proposition 8.18 implies that 
T^il^n) = for all n > 0. Proposition 8.17 and weak* convergence imply that 
ri{jl) for all i. Hence 
in 



= linXn^oo Tiijln) 

(see e.g. Theorem 2.1. 
that for all n > 1, 



Billingsley 



1999 



r^,(/i) = 0. The Portmanteau theorem 
applied to the closed set F^/^ implies 



K^i/n) > liminf/im(ri/„) 
> liminf/im,(ri/„) 



> £. 



^.15 



ll{nnTi/n) > £■ Remark 
)t'i w] 
1 and z>i(r+\ro) = 1. By Proposition [8l8 



implies there exist a > such that 



Therefore A(ro) 

/i = az^o + (1 ~ ct)^! where Dq, vi are G-invariant probability measures satisfying 



fiiyi) = for alH G {1, . . . , k}. 
Condition (a) implies that r^,{vQ) > 0, in which case = r^{fl) = ar^:{i>o) > which 
is a contradiction. ■ 



8.5. Proof of Theorem 8.10 , In the first part of the proof, we show that condition 



(a), (b) and 
Let $7 C be defined as in Lemma 



(a) implies the conclusion of Theorem 8.10 The second part is devoted to showing 
the equivalence of conditions 

^\S'o and 



-1/ 



8.20, Choose (uj'x') e Q x 



7 G TT^'^luj' ,x') C r+\ro. By definition of the set V, there exists a time T > 
such that ^^{uj',x') £ Q x V, for all t > T. Since // is a weak* limit point of 
the family {At{^'^ {lo' , x'))}t>o if and only if it is a weak* limit point of the family 
{At{uj', x')}t>o, we do not loss generality by considering {A((<I>^(a;', x'))}t- Since Q x 
V is compact, the set of all weak* limit points of the family of probability measures 
{At(<I>-^(a;',x'))}tgN is a non-empty subset of ^(17 x V). Let fi = liuik^oo-^tki^^^) 
be such a weak* limit point. Since r_|_ is positively 0-invariant and compact, pass- 
ing to a subsequence if necessary, there exists jl = limfe_!.oo Af^. (7) G 5'(r_|_). By 



Lemma 8.20, jj, G InvQ(0)(r+). Furthermore by Lemma 8.19 and continuity of ttq, 
^o(/^) ~ /"• Hence, Proposition 8.17, the continuity of the map ^, and property (ii) 



of Proposition 8.14[ imply the following equalities for all 



i: 



n{i^) 



Civ)dfii'n) 



ifc-1 



hm i Yl '^(®'(^)) 

fe— >oo Efc ^— ' 



s=0 



< ri{j). 



Hence, by Corollary 8.16 



rj(/2) < 0, for all i. 



Remark 8.15 implies there exists a > such that jl = ai/Q + (1 — a)i>i where Uj 
are invariant probability measure satisfying Dq{Tq) = 1 and t'i(r+\ro) = 1. By 
Proposition 8.18, rj(t'i) = for all i G {1, . . . , k}. Condition (a) implies r=i,(z>o) > 0. 
Therefore a must be zero, i.e. /i(r-|_\ro) = 1. Fix e > 0. By Lemma |8.21 there 
exists 77* > such that 
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which imphes 

fi{n X 5^) < e, Vr/ < rj* . 

Smce 7]* does not depend on /i, we have 



t—>oo 



for all x' G M" \5o and ui' G O, which implies that (ll) is stochastically persistent. 

Second, we show the equivalence of the conditions (a) and (b). We need the 
following version of the minimax theorem (see, e.g., Simmons [1998|): 



Theorem 8.22 (Minimax theorem). Let A, B he Hausdorff topological vector spaces 
and let L : A X B ^ M be a continuous bilinear function. Finally, let E and F be 
nonempty, convex, compact subsets of A andB, respectively. Then 



We have that 



min max L(a,b) = max min £ (a, b) . 
a£E beF beF a£E 



min max r j (/i) = min max > pirAjl) 



where the minimum is taken over fi € InV(Q(B)(ro) and the maximum over p S 
A := {p G M'p : X^jPi = 1}. Define A to be the dual space to the space 
bounded continuous functions from Fq to M and define B = M™. Let S = A, and 
D = InV(Q(0)(ro) C A which is nonempty, convex and compact by Proposition 
8.6, Let £ : j4 X i? — )■ M the bilinear function defined by L{jl,p) := Yli Pi''~iil^) ■ 
Proposition 8.17| implies that L is continuous. With these choices, the Minimax 
theorem implies that 



(30) 



min max r j (p) 

fj, i 



'E: 



max mm > piri{jl) 



where the minimum is taken over fl G InvQ)(0)(ro). By the ergodic decomposition 



theorem, the minimum of the right hand side of (30) is attained at an ergodic 



probability measure with support in Lq. Thus, the equivalence of the conditions is 
established. 

Third, we show the equivalence of condition (b) and (c). As a direct consequence 
of assertion (i) of Proposition 8.9, condition (c) implies ( b). T o prove the other 
direction, let fi C il be defined as in the proof of Lemma 8.20 Choose (a;',x') S 
Q.xSq. By the same arguments as above, there exist T > 0, 7 E ttq" (<I>"^(cj', x')) C Fq 
and pL G InV(Q(0)(ro) such that 



riifi) 



C{v)dfiiri) 






< nil)- 



Assertion (i) of Proposition 



8.9 



implies that rj(7) = rj(<l>-*^(a;', x')). Since <I>"^(a;',x') 
is on the same trajectory that (a;',x'), rj(/l) < rj(<l>^(u;', x')) = ri{io',x'). Writing 
/i as a convex combination of ergodic probability measures, condition (b) implies 
j:^Piri{uj',x')>0. U 
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